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ABSTRACT 


Current  ramp  designs  used  in  RORO  operations  have  been  determined  to  be 
structurally  inadequate  in  Sea-State  3.  The  main  reason  for  this  is  that  when  the  ramp  is 
loaded  with  two  vehicles,  the  relative  motions  between  the  ship  and  the  discharge  facility 
induce  stresses  above  yield.  The  objective  of  this  thesis  is  to  formulate  the  problem  and 
present  results  that  will  enable  the  design  of  proper  isolators.  The  problem  is  formulated 
in  terms  of  the  hydrodynamic  interaction  between  adjacent  bodies  along  with  structural 
coupling.  The  applicability  and  limitations  of  strip  theory  approaches  are  established 
through  comparisons  with  three  dimensional  hydrodynamic  analysis  data.  An  analytic 
model  of  the  ship-ramp-barge  system  is  developed  and  tested.  The  results  indicate  the 
validity  of  the  approach  and  establish  a  procedure  that  may  be  utilized  for  the  design  of 
passively  controlled  isolators. 
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I.  INTRODUCTION 


Over  the  years,  there  has  been  very  little  research  conducted  on  vessels  with 
objects  in  the  vicinity.  Of  the  research  conducted,  most  of  it  was  centered  on  the  analysis 
of  supply  ships  mooring  alongside  oil  platforms.  With  the  increased  emphasis  on  littoral 
warfare,  it  is  becoming  more  important  to  be  able  to  discharge  supplies  from  ships  to 
shore  without  a  readily  available  port  facility.  Therefore,  it  is  desired  to  develop  supply 
ships  with  the  ability  to  offload  cargo  to  a  roll-on,  roll-off  discharge  facility  (RRDF)  in 
the  open  ocean. 

To  accomplish  this,  a  ramp  is  mounted  from  the  ship  to  the  RRDF.  However,  due 
to  hydrodynamic  forces  on  the  ship  and  barges,  the  mounting  connections  are  exceeding 
their  design  limitations.  Therefore,  it  has  become  necessary  to  develop  a  method  to 
minimize  these  forces  on  the  pinned  connections.  To  accomplish  this,  the  hydrodynamic 
forces  on  the  ship  and  barge  must  be  accurately  modeled.  Then,  these  hydrodynamic 
forces  must  be  adapted  to  determine  the  magnitude  of  the  forces  that  act  on  the  pinned 
connection  so  that  a  passive  damping  system  may  be  produced  in  the  future. 
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II.  BACKGROUND 


A.  MOTIVATION 

In  recent  years,  the  United  States  military  has  taken  an  increasing  interest  in 
developing  the  ability  to  offload  equipment  and  machinery  from  roll-on,  roll-off  (RORO) 
vessels  in  the  open  seas.  Operational  requirements  for  equipment  transfer  are  set  at  sea 
state  three  while  maintaining  a  total  dynamic  acceleration  of  no  more  that  25%  of  the 
static  gravitational  constant.  In  order  to  accomplish  equipment  transfer,  a  connected  ramp 
system  has  been  established  between  the  RORO  and  a  transport  barge.  The  operational 
requirements  stated  above  should  hold  in  the  case  of  two  tanks  located  in  the  middle 
point  of  the  ramp.  This  corresponds  to  an  envisioned  worse  case  scenario  of  one  tank 
breaking  down  in  the  middle  of  the  ramp  and  the  second  tank  providing  assistance.  This 
coupled  problem  of  multiple  rigid  bodies  combines  elements  of  hydrodynamic 
interactions  and  structural  response.  The  philosophy  of  a  design  and  analysis  strategy  is 
how  to  ensure  proper  connectivity  between  the  different  bodies  such  that  cargo  motions 
and  ramp  stress  levels  are  minimized. 

B.  OBJECTIVE 

The  objective  of  this  thesis  is  to  formulate  and  present  preliminary  results  for  the 
problem  of  reducing  ramp  stress  levels  via  reductions  in  ramp  motion  inputs.  This  will  be 
done  by  coupling  the  hydrodynamic  analysis  of  the  adjacent  floating  bodies  to  a 
structural  model  of  the  ramp.  Passive/semi-active  isolator  design  will  then  be  used  as  the 
basis  for  the  synthesis  of  the  entire  model  in  the  frequency  domain. 
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C.  BACKGROUND 


A  great  deal  of  the  analysis  of  hydrodynamic  study  of  multiple  body  systems 
revolves  around  the  wave  effects  upon  vertical  cylinders  in  waves.  Most  of  this  research 
surrounds  the  desire  to  accurately  predict  the  affects  of  waves  on  different  parts  of  a 
structure  to  ensure  geometric  stability.  Much  of  this  research  is  motivated  to  produce 
continued  stability  in  offshore  oil  producing  stations  during  drilling  operations. 

However,  there  has  not  been  significant  analysis  performed  when  two  independent  bodies 
are  in  close  proximity  to  one  another.  [Ref.  9:  p.  725] 

Developing  the  equations  of  motion  of  a  solitary  vessel  in  various  incident  wave 
conditions  is  fairly  routine.  It  simply  involves  defining  the  added  mass  characteristics 
and  associated  stiffness  and  damping  matrices,  which  the  ship  will  possess  due  to  its 
design.  However,  the  effects  that  waves  will  have  on  ships  in  close  proximity  to  one 
another  are  very  different.  This  is  because  the  diffracted  and  radiated  waves  and  their 
respective  hydrodynamic  forces  will  affect  the  adjacent  vessels. 

The  equations  of  motion  of  vessels  in  close  proximity  to  one  another  are  derived 
in  much  the  same  manner  as  independent  vessels.  However,  when  multiple  vessels  are 
close  to  one  another,  additional  terms  must  be  incorporated  into  the  equations  of  motion 
to  ensure  validity.  The  importance  of  these  additional  terms  depends  on  relative  size  of 
each  vessel,  how  each  vessel  is  positioned  with  respect  to  the  other,  and  overall  distances 
between  them.  [Ref.  9]  Taking  these  criteria  into  account,  an  additional  added  mass  term 
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and  an  additional  damping  term  must  be  utilized  to  accurately  predict  the  hydrodynamic 
response  of  each  vessel. 

Research  has  indicated  that  there  are  multiple  techniques  that  could,  in  principle, 
predict  the  numerical  values  of  the  additional  terms  necessary  to  accurately  predict  the 
hydrodynamic  motion  of  a  vessel.  One  of  these  techniques  is  the  use  of  what  is  known  as 
the  macroelement  technique.  This  method  approximates  each  of  the  vessel’s  cross- 
sectional  areas  with  a  step  curve.  The  flow  field  around  each  area  is  then  subdivided  into 
macroelements  and  approximated.  This  technique  uses  the  method  of  separation  of 
variables  and  Laplace’s  differential  equations  for  the  flow  field.  Once  the  velocity 
potentials  for  the  diffraction  and  radiation  features  for  each  element  are  established,  a 
Fourier  series  is  developed  that  meets  the  free  surface  and  ship  body  boundary 
conditions.  To  obtain  the  remaining  coefficients,  the  Galerkin  method  is  utilized.  This 
method  is  effective  as  it  does  not  produce  many  numerical  irregularities  and  saves  in 
computer  processing  time.  [Ref.  1]  Due  to  its  approximations,  this  method  may  be 
inaccurate  for  certain  cases  and  it  suffers  from  the  basic  limitations  of  a  two-dimensional 
analysis,  most  notably  end  effects  and  the  existence  of  standing  waves. 

Another  technique  that  is  widely  studied  is  that  of  finite  element  approximations. 
In  this  approach,  the  area  around  the  vessel  is  modeled  using  boundary  dampers/infinite 
elements  or  eigenfunction  expansion/boundary  integral  techniques.  [Ref  5:  p  358]  Using 
standard  assumptions  that  wave  and  body  motions  are  small,  the  complex  velocity 
potential  is  obtained  for  the  diffracted  and  radiated  wave  potentials.  Using  data  regarding 
complex  velocity  potentials;  information  about  linearized  free  surface,  bottom,  radiation. 
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and  body  boundary  conditions;  and  hydrodynamic  forces  and  coefficients,  a  final 
equation  is  developed.  This  equation  integrates  and  sums  the  added  mass,  damping,  and 
stiffness  matrices  of  each  vessel  and  its  effects  on  the  other  vessel.  Furthermore,  it  solves 
for  the  forces  on  each  vessel  due  to  waves  and  any  associated  hydrodynamic  effects  from 
the  other  vessel.  This  method  is  effective  for  various  wave  conditions  as  long  as  the 
problem  assumptions  are  met.  However,  proper  nodal  placement  and  meshing  is 
essential  for  accurate  results.  Due  to  the  requirement  for  very  fine  meshing  in  some 
problems,  computer  resources  may  restrict  the  determination  of  accurate  results.  [Ref  5] 

Techniques  based  on  strip  theory  formulations  are  easier  to  solve  than  the 
previous  methods  while  producing  an  accurate  wave  force  and  linked  body  prediction.  In 
this  approach,  there  are  two  floating  vessels  attached  alongside  one  another.  The 
equations  of  motion  are  generated  using  known  data  for  each  vessel  and  coupled  together. 
When  coupled,  each  equation  of  motion  from  one  vessel  effects  that  of  the  other.  Vessel 
interaction  is  approximated  by  assuming  that  as  a  wave  strikes  a  vessel,  the  wave  is 
reflected  an  infinite  number  of  times  between  vessels.  This  produces  two  additional 
equations  for  the  wave  interaction  forces  on  each  vessel.  To  determine  the  wave 
interaction  forces,  strip  theory,  which  assumes  that  changes  in  the  longitudinal  gradients 
of  the  flow  field  around  a  vessel  are  much  smaller  than  corresponding  changes  in  the 
transverse  direction,  is  utilized.  The  wave  forces  are  a  combination  of  the  Froude-Krylov 
force  and  the  diffraction  wave  force.  Additional  equations  for  added  mass,  damping,  and 
stiffness  for  each  vessel  are  needed  to  represent  the  interaction  effects.  Taking  the 
equations  generated  by  wave  interaction  between  the  vessels  and  those  equations  for 
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added  mass,  damping,  and  stiffness  and  placing  them  back  into  the  original  equation  of 
motions  for  each  vessel,  a  new  equation  of  motion  is  developed  which  will  depict  the 
motion  of  the  vessel  to  various  wave  loading.  This  technique  has  shown  to  be  effective 
in  all  sea  directions  as  long  as  the  assumptions  of  linearized  motions  and  slender  bodies 
are  accurate.  [Ref  2] 

Another  related  technique  that  is  commonly  used  to  solve  for  the  additional  terms 
is  that  of  the  boundary  integral  method.  This  method  first  divides  the  vessels  into 
sections.  The  radiation  velocity  potential  function  is  determined  using  Green’s  theorem 
and  the  Fredholm  integral  equation  of  the  second  kind,  which  is  then  applied  to  each  of 
the  divided  sections.  Once  the  radiation  velocity  potential  function  is  developed,  an 
interpolating  shape  function  can  be  formed.  Finally,  a  set  of  linear  equations  are 
developed  which  produce  a  set  of  algebraic  equations  for  each  of  the  radiation  velocity 
potential  equations  throughout  the  vessels.  This  method  produced  accurate  results  in  all 
types  of  wave  loading  directions  at  varying  frequencies.  [Ref  3] 

D.  CONTENTS 

The  boundary  integral  technique  is  the  method  that  most  researches  utilize  when 
obtaining  the  numerical  solution  to  the  added  mass,  damping,  and  stiffness  terms  in  the 
equation  of  motion  of  a  vessel.  For  this  reason,  it  was  adopted  as  the  method  of  choice  for 
this  study.  Chapter  III  presents  the  theoretical  foundation  for  the  hydrodynamic  equations 
of  motion.  This  is  followed  by  the  development  of  the  algebraic  constraint  equations  that 
were  incorporated  in  the  original  differential  equations  of  motion  in  order  to  model  the 
structural  coupling  between  the  bodies.  Modeling  of  the  connecting  ramp  is  also 
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included  in  Chapter  III,  followed  by  the  results  of  this  work  in  Chapter  IV.  Finally,  the 
conclusions  from  this  research  and  some  recommendations  for  follow-up  work  are 
summarized  in  Chapters  V  and  VI. 
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m.  PROBLEM  FORMULATION 


A.  INTRODUCTION 

When  examining  the  connection  between  the  ship  and  the  barge,  the  points  where 
the  ramp  attached  to  the  barge  and  the  ship  exhibit  induced  forces  on  the  pinned 
connection.  The  equations  to  develop  these  forces  are  formulated  by  considering  heave, 
pitch,  and  roll  forces  upon  the  ship  and  the  barge. 

B.  HYDRODYNAMIC  EQUATIONS  OF  MOTION 

To  obtain  the  force  equations,  it  was  necessary  to  define  a  standard  coordinate 
system.  The  system  used  was  the  same  one  utilized  in  the  reference  14,  Principles  of 
Naval  Architecture,  in  which  r|i=surge;  r|2=sway;  r)3=heave;  rj4=roIl;  r)5=pitch;  and 
r)6=yaw.  The  three  principle  motions  that  affect  the  problem  were  heave,  roll,  and  pitch. 
Relative  surge,  sway,  and  yaw  motions  were  considered  negligible  and  were  not  analyzed 
further. 


The  other  assumption  that  was  made  at  this  point  was  that  the  only  appreciable 
affect  upon  ship  movement  would  be  in  the  vertical  direction.  Therefore,  the  equation  for 
absolute  vertical  motion  of  a  vessel  at  a  given  point  was  found  to  be 

Is  =  +  n4y-xTi5  (l) 

Where  yand  x  is  the  motion  at  any  point  on  the  body  in  the  x  and  y  planes.  £3  is  the 
complex  local  amplitude  at  x,  and  y  in  heave. 

After  defining  the  governing  equation  that  would  control  the  motion  of  each 
vessel,  a  ffee-body  diagram  was  produced  of  the  ramp.  This  diagram  assumed  that  there 
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are  three  spring  connections  mounting  the  ramp  to  the  ship  and  barge.  Each  connection 
produces  a  displacement  at  the  pinned  joint  defined  as  £1,  £2,  and  £3.  The  first  £ 
represents  the  connection  for  the  ship  and  the  last  two  £’s  are  for  the  barge  connections. 
The  next  assumption  was  that  each  spring  would  exert  a  force  at  the  connection  as  a 
force-couple  system.  This  force-couple  system  represents  a  single  force  at  each  pinned 
location.  Therefore,  the  equation  given  for  the  force  at  each  location  is  defined  as 

F,  =  G,(, 


(2) 


where  the  subscripts  i  and  j  reference  the  particular  force-couple  system  at  a  given 
location.  G  is  the  matrix  of  influence  coefficients,  which  is  a  transfer  function.  From  (2), 
the  following  matrix  was  obtained: 


'K 

Gn 

Gn 

Gn 

h' 

~ 

G2l 

g22 

G2i 

£ 2 

A. 

.^3. 

G32 

Gn_ 

A. 

Using  this  process,  it  was  then  necessary  to  find  the  equations  of  motion  for 
heave,  roll,  and  pitch.  This  was  done  using  (4)  where  A,  B,  and  C  are  the  added  mass, 
damping,  and  stiffness  matrices  for  the  ship  and  barge.  The  j  terms  are  for  the  6 
directions  of  motion  in  which  only  heave,  roll,  and  pitch  were  considered.  Fj1  and  Fj°  are 
the  complex  amplitude  of  the  exciting  force  components  due  to  incident  and  diffracted 
waves,  and  F/  are  the  connection  forces  at  the  ramp,  to  be  determined. 

t  [-  (a, + 4. )+  +  c,  tT-  f' + f;  *f;  ( 4 } 

Solving  these  equations  for  the  three  relevant  motions  to  this  problem,  the  final  equations 
of  motion  were  found  to  be: 
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AssTj  ^  +  Az$rjs  -  F2  in  Heave  '  * 

A  44  tj  4  =  F4  in  Roll  (  6  ) 

As3Tj3  +  As5ijs  =  Fs  in  Pitch  (  7  ) 

The  various  terms  in  these  equations  are  defined  as: 

77  3  (f)  =  The“°‘'  =>  Heave  (  8  ) 

rj4(t)=T]4e,a'1  =>  Roll  (9) 

75(f)=  775e'“'' =>  Pitch  (10) 


Each  of  the  internal  terms  in  the  equations  of  motion  are  all  matrices.  These  matrices  are 
characterized  as  A,  B  and  C  and  represent  the  added  mass  of  the  vessel,  the  damping 
coefficient,  and  the  stiffness.  This  implies  that  each  vessel’s  hydrodynamic  properties  are 
modeled  as  a  spring-mass-damper  system.  The  first  number  in  the  following  equation’s 
subscripts  refers  to  the  direction  (heave,  pitch,  or  roll)  and  the  second  number  is  with 
respect  to  a  particular  direction.  For  example,  A35  in  (12)  would  be  the  heave  added  mass 
due  to  pitch. 


A33  --^e2(A+A33)+ja)eB}3  +c33 

(ii) 

As  =-tfj2( Axe  +A3^+ia)iiB}5+C3S 

(12) 

Ass  =-0)c2(I55  +A55)+io)eB55  +C55 

(13) 

As  =-#>2(AXc  +^3)+/fi^j^3+Q3 

(14) 

Am — — ty  2(/44  +4^) +^,£44  +C44 

(15) 

The  terms  in  the  added  mass  equations  are  the  specific  equations  of  mass,  damping,  and 
stiffness  for  the  particular  vessel.  Due  to  inherent  symmetrical  relations  within  the 
vessel,  the  following  additional  assumptions  are  utilized: 
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(16) 

(17) 

(18) 


C 

B 

A 


53 

35 

35 


A  53 


The  force  terms  in  the  equations  of  motion  are  simplified  to  make  individual  force 
calculations  easier.  This  leads  to  the  following  equations  that  produce  the  resultant 


heave,  roll,  and  pitch  forces: 


F  3  =  F3‘  +  F3d 

+  F3r 

(19) 

F  4  =  F  /  +  Fad 

+  F/ 

(20) 

F  5  =  F5'  +  F  SD 

+  F  $r 

(21) 

The  F3r,  F/,  and  Fsr  terms  are  the  heave,  roll,  and  pitch  forces  that  act  on  the  ship  or  barge 
from  the  motion  of  the  ramp. 

The  equations  of  motion  in  (5)  through  (7)  are  body  independent.  Consequently, 
they  are  used  to  describe  the  motions  of  the  ship  and  barge.  Therefore,  a  subscript  ‘S’  or 
‘B’  is  used  to  distinguish  the  motions  of  each  vessel.  Furthermore,  equations  (5)  through 
(7)  fail  to  recognize  that  as  two  objects  are  in  close  proximity  to  one  another,  there  will 
exist  a  radiative  effect  between  the  two  vessels.  This  results  from  waves  striking  one 
vessel,  bouncing  off  the  hull  and  then  striking  the  vessel  alongside.  These  reflected 
waves  effect  heave,  pitch  and  roll  of  the  adjacent  vessel.  Combining  the  vessel 
distinguishing  subscripts  and  the  radiative  matrices  yielded  the  following  equations. 

Heave,  ship  =>  A33,stj3S  +  A35,sT)iS  +  Bn,BTji  B  +  Bis,bT]}  B  +  Bu.bT]4B  =  F1S 
Pitch,  ship  =>  As3.sT]3S  +  Ass,stjss  +  B$3,btj3B  +  Bss,b1]5B  +  Bsa.bT]a  b  =  Fss 
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(22) 

(23) 


Roll,  ship  =>  Aaa,s  tj4S  +  Baxb7]3B  +  Bas,bT]sb  +  Baa,bti4B  =  ^4,5 

Heave,  barge  =>  Ai3,btj3B  +  A3s,btisb  +Bn,sT}iS  +  B35.sT]5S  +Bm,sT]4S  =  F3B 


Pitch,  barge  =>  As3,btj3B  +  As5,btj}B  +  Bs3.stj3S  +  Bss,sr/ss  +  Bs4,sT]4S  -  Fi'S 


Roll,  barge  A44,bT]4B  +  BA3.srj3S  +  B4s.stj5S  +  Baa,stj4S  =  F4B 


(24) 

(25) 

(26) 
(27) 


The  B  terms  are  the  radiative  matrices  between  the  barge  and  ship  while  the  subscript  S 
and  B  designate  the  ship  or  barge.  The  radiative  matrices  are  defined  with  the  following 
equations: 


8  33,5 

”  ®  ^33,5  +  icoE  b 

(28) 

B  34  ,5 

—  —co  D  34  B  +  i  (oE34  b 

(29) 

-835,5 

—  —co  D35  b  +  icoE3S  B 

(30) 

B  43,5 

—  —co  D43  B  +  icoE  43 B 

(31) 

-844,5 

=  —co  D 44  B  +  icoE 44  B 

(32) 

845,5 

=  -^^45,5  +  i0)E45,B 

(33) 

8  53  ,5 

=  -°)2D52  B  +  icoE53  B 

(34) 

854,5 

—  —co  D5 4  B  +  icoES4  B 

(35) 

8  55,5 

—  —CD  D 55  B  +  icoE55  B 

(36) 

D  and  E  are  the  added  mass  and  damping  influence  coefficients  due  to  the  radiated  waves 
from  one  vessel  to  another.  Similar  equations  for  the  radiative  matrices  are  obtained  for 
the  ship.  Due  to  symmetry  Ejk(1)=Ek/2)  and  Djk(1)=Dk/2).  An  example  of  one  of  these 
symmetric  equations  is  the  following: 

-^35,5  =  D53S  (37) 
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Using  a  free-body  diagram,  (5)  through  (7)  were  expanded  allowing  the 
hydrostatic  and  ramp  forces  on  each  vessel  to  be  separated.  The  terms  Fir,  F2r,  and  F3f 
were  defined  as  fi,  f2,  and  f3  respectfully.  The  following  revised  equations  of  motion  are: 


Heave,  ship=>  A33,s7j3S  +  A3s,stj5S  +  B33,b7]3B  +  B3a,bT]4B  +  £35,5775#  =  F3/+fl  (38) 

Pitch, ship:=>  A$3.s7]3S  +As5,st}5S  +B$3,btj3B  +B5a,btj4B  +BssjT]5j  -  F5Sh  -fxxx  (  39  ) 

Roll, ship  Aaa,s  TJas  +  Ba3,b  7J3B  +  Baa,b  7]a  b  +  £45,5  7]SJ  -  FJ*  +/>'!  ( 40  ) 

Heave,  barge  =>  A33,btj3B  +  A3s,btjsb  +  B33,stj3S  +  B3a,stj4S  +  B3s,stjss  =  Fi.BH  +fi+fi  (41) 
Roll, barged  Aaa.bt, 74  8  +Bn,sr)is  +  Bu,s7jiS  +B45.stj5S  =  FABH  +f2y2  +f2y2  (  42  ) 

Pitch, barge  =>  As3,b7]3B  +  Ass,Brj5B  +  Bs3,stj3S  +  Bsa,s7]4S  *4*  Bs$,stjss  =  f$.bH  (43) 


This  allowed  for  the  effects  on  the  ship  and  barge  at  the  pinned  connections  due  to  the 
direct  forces  from  the  ramp  as  well  as  moment  forces  due  to  ramp  movement.  The 
superscript  H  is  to  show  the  force  acting  on  the  hull  of  the  respective  vessel. 

Once  the  above  hydrodynamic  equations  of  motion  are  solved,  the  values  for  fj, 
f2,  and  fj  can  be  established.  However,  all  of  the  q  terms  in  the  above  equations  are 
unknown  and  all  the  f  terms  are  also  unknown.  Therefore,  the  values  for  the  q  terms  are 
determined  by  defining  additional  equations,  q  is  thus  defined  as  a  function  of  fi,  f2,  and 
fi  for  the  ship  and  barge  as: 
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*7  3  ,S  —  f*  31  f  l  +  M  30 

(44) 

*7  4,S  =-  t1  41  f  \  +  f*  40 

(45) 

*7  5,S  =  L1  51  f  \  +  M  50 

(46) 

V  3  ,B  =  V  32  f  2  ^  33  f  3  ^30 

(47) 

*7  4,S  ~~  ^  42  f  2  ■*"  ^  43  f  Z  ***  ^  40 

(48) 

*7  5,B  =  °.S2  f 2  +  °  53  f  %  +  V  50 

(49) 

From  these  equations,  it  is  possible  to  solve  the  six  linear  equations  for  each  respective  p 
and  v  in  terms  of  various  known  added  mass  matrices,  known  physical  positions  on  the 
vessel,  and  the  forces  acting  on  the  pinned  connections  as  well  as  the  hydrodynamic  wave 
forces. 

Once  each  p  and  v  term  is  solved,  every  r\  is  considered  a  known  value  in  terms 
of  fi,  f2,  and  f3.  Each  r|  term  has  a  contribution  to  the  total  force  acting  where  the  ramp 
attaches  to  the  barge  or  the  ship.  Therefore,  the  total  force  acting  on  a  pinned  joint 
includes  forces  from  heave,  roll,  and  pitch.  This  produces  three  combined  displacement 
equations  for  each  connection  that  are: 

#1  3,s  +  V*,sy\  -  x\Vs,s  =  kuf\  +  *,0  •  (  50  ) 

%2  =  V  3, B  V  4,By  2  ~  X2  V  i,B  =  k2lfl  kH  fi  ^20  (  51  ) 

~  V  3,B  *7  4,fi^3  —  X2r!  5  ,B  ~  k32  f 2  k33  fi  k30  (  ) 

It  was  stated  previously  that  and  £3  are  the  combined  displacements  that  act 

at  each  of  the  pinned  connections.  Now  the  equations  are  defined  by  relating  file  forces 
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acting  on  the  pinned  connections  from  the  above  k  terms  using  the  known  formulas  for  r\ 
and  equating  portions  of  the  £  equations  that  correspond  to  the  respective  x  or  y  terms. 
Thus,  the  k  terms  become: 


~  /*3i 

klo=M3o+M^y,-Msoxl 

(53) 

k22  =^32  +U«J'2“y52*2 

k23=u33  +  o^y2-u53x2 

^20  "  ^30  +  ^40^2  ^50*2 

(54) 

^32  =  ^32  ^42^3  ~  ^52*3 

ki3  =  t>33  +  vAJy3  —  oS3x 3 

^30  =  ^30  +  ^40^3  ”  ^50*3 

(55) 

These  k  terms  are  then  set  equal  to  the  transfer  functions  that  relate  the  forces  at  one 
pinned  connection  to  the  other  forces  from  each  pinned  connection,  and  how  the  various 
forces  interact  through  the  ramp  system.  This  leads  to  the  equations: 

^10  ^11  f\  ~  G II  f\  ■*"  ^12  fl  ^13  fs  (  56  ) 

^20  ^22  f  2  ^ 23  f  3  ~  G  2)  f\  +  G  22  f  2  +  G  23  f  3  (  57  ) 

^30  ■*"  ^32  f  2  "*■  ^33  f  2  =  G  31  f\  +  G  32  f  2  +  G  33  f  3  (  ^  ) 

Rearranging  equations  (56)  through  (58)  and  placing  into  matrix  for  yields: 


^•u  Glt 

“G12 

~G13 

X 

“*io 

-G2l 

^22  ^22 

^23  _^23 

/* 

= 

"^20 

.  -C3, 

^32  “^32 

*33  "^3. 

A 

_”^30_ 

This  first  matrix  is  labeled  as  matrix  [Gk].  Solving  for  the  forces  yields: 


f  2 

=  [Gj" 

*  10 
“  *  20 

f  3 

I 

o 

1 

(60) 


The  three  forces  are  easily  solved  using  Matlab  or  any  other  engineering  computational 
program  given  the  numerical  values  of  the  added  mass,  damping,  and  stiffness  matrices. 
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C.  RAMP  EQUATION  OF  MOTION 


1 .  Introduction 

To  accurately  model  the  ship  barge  combination,  it  is  important  to  incorporate  the 
ramp  connecting  the  two  vessels  because  there  will  be  interaction  effects  between  the  two 
vessels  transmitted  through  the  ramp.  The  ramp  is  modeled  here  as  a  rigid  body, 
although  the  procedure  is  general  enough  so  that  it  can  be  generalized  to  include 
flexibility  effects. 

2.  Derivation  of  the  Equation  of  Motion 

To  obtain  the  equations  of  motion  for  the  ramp,  it  was  first  necessary  to  examine 
the  ramp  as  a  separate  entity.  First,  the  ramp  is  a  three-degree  of  freedom  system  in  the 
heave,  pitch,  and  roll  directions.  The  ramp  is  assumed  to  be  connected  to  the  ship  and 
barge  at  three  locations.  See  appendix  A  for  figure  of  LMSR/barge  connection.  Each  of 
these  locations  has  a  spring  connection  and  the  overall  ramp  has  a  mass  of  M.  It  was 
further  assumed  that  the  ramp  is  a  homogenous  device  with  uniform  properties 
throughout.  Using  these  assumptions,  the  equations  of  motion  were  derived  using  the 
method  of  summing  kinetic  and  potential  energies  of  the  ramp  and  then  using  Lagrange’s 
Equations  to  derive  the  equations  of  motion. 

First,  it  is  necessary  to  define  Lagrange’s  equations.  The  Lagrangian  is  defined  as 

L  =  T  - U  (61) 

where  L  is  the  Lagrangian,  T  is  kinetic  energy,  and  U  is  the  potential  energy  of  the  ramp. 
The  Lagrangian  is  then  used  to  find  the  value  of  the  forces  at  each  degree  of  freedom. 
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The  equation  for  this  is  written  below  as  equation  (62).  This  equation  takes  the  partial 
derivative  of  the  Lagrangian  with  respect  to  the  velocity  of  one  of  the  degrees  of  freedom. 
The  partial  derivative  of  the  Lagrangian  with  respect  to  position  is  then  taken  and 
subtracted  from  the  time  derivative  of  the  first  value.  This  produces  the  force  at  one  of 
the  degrees  of  freedom. 


d_ 

dt 


f  \ 
8L 

d  <h  J 


8L_ 

dq 


=  F, 


(62) 


After  the  Lagrangian  is  defined,  the  kinetic  and  potential  energies  for  the  ramp  are 
found.  Looking  at  the  ramp  from  a  view  approximately  45  degrees  above  the  horizontal, 
it  was  assumed  that  there  are  three  pinned  connections.  These  are  at  the  two  comers  that 
would  attach  the  ramp  to  the  barge  and  the  third  connection  is  from  a  kingpin  in  the 
center  of  the  ramp  that  attaches  the  ramp  to  the  ship.  These  connections  produce  a  force 
upward  due  to  the  associated  springs.  Visualizing  the  ramp  in  some  deformed  state,  it  is 
possible  to  define  the  change  in  position  due  to  heave,  pitch,  and  roll.  The  heave  is 
simply  the  distance  y  upward  from  the  original  axis.  The  pitch  is  defined  as  a  positive 
upward  angle  from  the  new  horizontal  axis  on  the  ship  side  of  the  ramp.  Roll  is  the 
rotation  about  the  x-axis  from  the  new  pitch  angled  ramp.  Kinetic  energy  is  defined  as 


2  2  2 
T  =  .5 M  y  +  .5M  6  +  .5M  w 

In  this  equation,  M  is  the  mass  of  the  ramp  and  y,  0,  and  T'  are  the  three  individual 
displacements  which  act  on  the  ramp  throughout  the  ramp  motion. 
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Potential  energy  is  defined  as  the  stored  energy  of  the  springs  throughout  the 
system.  Damping  of  the  ramp  was  not  included  in  the  determination  of  the  potential 
energy  as  proportional  damping  was  considered  to  be  in  effect  Therefore,  the  value  of 
the  damping  coefficients  will  be  proportional  to  the  spring  coefficients  by  some  set 
amount.  This  allows  easier  calculation  of  the  potential  energy.  In  the  potential  energy 
equation,  U  is  the  potential  energy,  k  is  the  spring  coefficient,  A  is  the  length  of  the  ramp, 
and  B  is  the  width  of  the  ramp.  Y,  0,  and  'P  are  the  same  as  previously  defined. 

U  =  .5&,^y  +  ysin©j  +.5^2^-ysin©  +  ysin  'pj  +  .5£3^y-ysin©-ysin'pj  (64) 


Once  the  kinetic  and  potential  energies  are  defined,  the  determination  of  the 
equations  of  motion  of  the  system  is  fairly  straightforward.  Set  up  the  Lagrangian,  take 
derivatives,  and  equate  terms  and  place  in  matrix  form.  The  force  equations  were  found 
to  be 


A.^.B. 


M  y-k\  y  H - sin  ©  -  k2\  y - sin  0  H - sin  'F  -  k3 \  y - sin  © - sin  ¥  =  F 


A  .  ^  B  . 


(65) 


M&-— cos©&,l  ,y+— sin©!-— cos©&2l  y~~ sin©+— sin'F  cos©^!  j-— sin©-— sin'Fj  =  Fe  (66  ) 


A  .  ^  B  . 


cos'¥k2  y - sin  ©  +  — sin  ¥ - cos  Wk3  y - sin© - sin  T 


(67) 


It  should  be  noted  that  there  is  not  a  term  for  T  in  the  ki  part  of  the  equation.  This  is 
because  as  the  ramp  rolls,  the  changes  in  the  'P  angle  for  the  spring  on  the  centerline  is 
zero.  This  also  causes  the  Fy  force  to  not  have  a  ki  term  as  the  derivative  does  not  have  a 


'P  term  in  it. 
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(68) 


Equation  (68)  above  shows  the  end  result  matrix  for  the  three  directional  forces  on  the 
ramp.  This  matrix  is  produced  after  deriving  the  equations  of  motion  for  the  ramp, 
adding  the  proportional  damping  values  and  linearizing  the  all  three  equations  using  the 
small  angle  approximation.  This  allowed  sin0  and  sin¥  to  be  approximated  as  simply  0 
and  ¥  and  cos0  and  cos¥  to  be  approximated  as  1. 


Once  the  final  matrix  for  the  ramp  equations  of  motion  was  derived,  these  forces 
representing  the  position  of  the  ramp  in  terms  of  y,  0,  and  ¥  had  to  be  correlated  to  be  in 
terms  of  Fj,  F2,  and  F3,  the  forces  at  each  of  the  connecting  locations  and  their  associated 
changes  in  absolute  vertical  motion.  First,  the  forces  on  the  ramp  were  written  in  terms 
of  Fi,  F2,  and  F3.  This  produced  the  following  matrix  where  all  symbols  are  as  previously 
stated.  f  1 


(69) 


The  next  step  was  to  state  the  changes  in  position  of  the  ramp  connecting  locations  in 
terms  of  y,  0,  and  ¥. 


0 

_ B_ 

2  2 
-A  -  B 


2 

-  A 


(70) 


In  this  matrix,  xi,  X2,  and  X3  are  the  vertical  displacements  for  each  of  the  pinned 
connections.  The  final  relation  must  be  to  have  Fy,  Fe,  and  Fy  in  terms  of  only  the 
displacements.  The  resulting  equation  is  as  follows: 


A 


=  [-fi)2[M]+  ;©[c]+[a:]|  © 
20 


(71) 


In  this  matrix,  the  M,  C,  and  K  matrices  are  the  mass,  damping,  and  spring  stiffness 
matrices  defined  earlier.  Frequency  is  defined  as  to  and  is  the  frequency  of  wave  impact 
from  the  hydrodynamic  information. 


Once  the  relations  are  made  between  the  pinned  connections  and  the  equations  of 
motion,  a  relation  can  be  established  to  relate  the  displacements  of  the  connections  to  the 


ramp  equations  of  m 


‘  xl ' 

1  —  0 

2 

1  1 

1  1 

*2 

— 

1  ^  5 

2  2 

i  -B 

2  2 

[F»,a,  r 

A  - A 

2  2 

0  — 

L  2 

1  1 

to  X 

i_ 

Fi 

Fz 

f3 


(72) 


Known  values  are  the  length  and  width  of  the  ramp,  the  M,  C,  and  K  matrices,  and  the 
hydrodynamic  forces  on  the  pinned  connection.  Placing  all  this  information  together 
results  in  the  final  matrix  describing  the  change  in  vertical  displacement  of  the  pinned 
connection. 
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IV.  RESULTS 


A.  LARGE  MEDIUM-SPEED  ROLL-ON  ROLL-OFF  VESSEL  ANALYSIS 


1.  Introduction 

The  comparison  between  the  strip  theory  and  full  three-dimensional  (WAMIT) 
program  solvers  are  presented  in  Figures  1  through  48.  Each  of  the  values  in  the  figures 
has  been  non-dimensionalized  in  a  consistent  manner  displayed  in  Table  1.  Furthermore, 
all  pertinent  ship  parameters  are  shown  in  Table  3.  All  symbols  used  are  identified  in 
Table  2. 


Ship/Barge  Coefficient 

WAMIT 

Non-dimensionalization 

A33,  Heave  Added  Mass 

A^/pL? 

B33,  Heave  Damping 

B33/copLJ 

A55,  Pitch  Added  Mass 

A^L5 

B55,  Pitch  Damping 

B55/ CD  pi/* 

A3 5,  Heave  to  Pitch  Added  Mass 

A35/pL4 

B35,  Heave  to  Pitch  Damping 

B35/(opL4 

A53,  Pitch  to  Heave  Added  Mass 

A53/pL4 

B53,  Pitch  to  Heave  Damping 

Bss/opL* 

C33,  Heave  Stiffness 

C33*pgAw/pgL‘: 

F3,  Heave  Exciting  Force 

Fj/pgAL4 

F5,  Pitch  Exciting  Moment 

Fj/pgALJ 

Phase  Angle 

None 

Table  1.  Table  of  Various  Non-Dimensionalizations 


2.  Data  Analysis 

Figures  1  through  28  display  the  added  mass  and  damping  results  for  the  large 
medium  speed  roll-on  roll-off  vessel  (LMSR)  data.  The  first  8  figures  portray  data 
obtained  from  the  strip  theory  solver  and  WAMIT  when  over  2000  panels  were  used  in 
the  3-D  analysis.  The  next  8  figures  show  the  same  added  mass  and  damping  curves 
comparing  strip  theory  and  WAMIT  but  in  these  curves,  there  were  only  about  1000 
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panels  for  the  ship  and  barges.  Figures  21  through  28  show  the  effects  due  to  wave 
radiation  effects  from  the  barge  on  the  LMSR.  Overall,  the  results  show  that  relatively 
good  agreement  is  obtained  between  the  two  computer  solving  techniques,  especially  in 
the  heave  mode.  However,  oscillatory  behavior  is  observed  for  the  WAMIT  program  at 
high  frequencies  (low  periods),  especially  when  fewer  panels  are  utilized.  This  is  due  to 
discretization  effects  in  the  solver.  To  overcome  these  problems,  a  relatively  high 
number  of  panels  are  required  to  accurately  represent  the  hydrodynamic  interactions 
between  the  LMSR  and  the  small  Roll-on,  Roll-off  Discharge  Facility  (RRDF).  The  strip 
theory  solver  does  not  take  such  interactions  into  consideration.  Therefore,  it  is  not 
susceptible  to  high  frequency  “chattering”,  unlike  the  3-D  solver. 


SYMBOL 

DEFINITION 

M 

Mass  of  the  ship  or  barge  in  pounds. 

L 

Length  of  the  ship  in  feet. 

g 

Force  due  to  gravity  in  pounds-force  per  second-squared. 

Aw 

Waterplane  area  produced  by  the  strip  theory  solver  in  square  feet. 

A 

Wave  amplitude. 

G) 

Frequency  of  wave  impact  upon  the  ship  or  barge  in  Hertz. 

V 

Ship’s  displacement  in  cubic  feet. 

p 

Density  of  Seawater  in  pounds-mass  per  cubic  feet.  Considered  to  be  1.9904  lb/ft3. 

Table  2.  List  of  Applicable  Symbols  and  Their  Definitions 


As  mentioned,  the  3-D  solver  can  overcome  the  high  frequency  oscillations  by 
utilizing  a  more  refined  panel  scheme.  However,  the  very  large  number  of  panels  that  is 
required  by  the  3-D  solver  incurs  a  significant  cost  in  both  computational  time  (more  than 
16  hours  for  the  2000  panel  case  on  a  high-end  CPU)  as  well  as  the  time  needed  to 
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prepare  the  data  file.  By  contrast,  the  strip  theory  solver  results  were  obtained  after  only 
a  few  seconds  of  execution  time  on  an  older  generation  Pentium  CPU.  Of  course,  the 
limitations  of  the  strip  theory  solver  are  that  the  interaction  effects  are  not  capable  of 
being  modeled.  By  comparing  the  strip  theory  and  WAMIT  solvers,  it  was  clear  that 
modeling  of  the  added  mass  and  damping  characteristics  of  the  main  ship  could  be 
achieved  both  accurately  and  efficiently  through  a  combination  of  the  two-dimensional 
and  three-dimensional  hydrodynamic  solvers. 

The  results  for  the  hydrodynamic  wave  excitation  forces  and  phase  angles  for  the 
LMSR  are  presented  in  Figures  25  through  60.  These  wave  forces  range  from  zero 
degrees  (head  seas)  through  315  degrees  at  45-degree  intervals.  This  provides  complete 
360-degree  coverage  for  wave  force  analysis  for  the  LMSR.  In  these  plots,  overall 
agreement  between  the  strip  theory  solver  and  WAMIT  is  outstanding.  As  will  be  noted 
in  more  detail  when  the  barge  data  is  analyzed,  there  are  some  slight  differences  between 
the  port  and  starboard  wave  directions  in  the  3-D  results.  These  deviations  are  due  to 
interactions  between  the  barge  and  the  LMSR. 


Dimension 

LMSR 

Barge 

Length  (in  Feet) 

894.53 

79.2 

Beam  (in  Feet) 

105.8 

72.6 

Depth  (in  Feet) 

None  Given 

4.7 

Draft  (in  Feet) 

27.5(FWD),  32.27(AFT) 

1.3 

Displacement 

45,211  Long  Tons 

402,783  Pounds 

Ramp  Location  (Feet  From  FP) 

X=301.2,  Y=71.3 

39.5  (From  Hinged  Side  Edge) 
8.8(From  Centerline  of  Causeway 
to  Port) 

Table  3.  Applicable  Ship,  Barge,  and  Ramp  Dimensions 
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3.  Figures 


x  to’3  Heave  Added  Mass  -  LMSR 


Figure  1.  Plot  of  Strip  Theory  and  Wamit 
Program  Comparisons  for  Heave  Added  Mass 
for  the  LMSR. 


x  io"3  Heave  Damping  -  LMSR 


Figure  2.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Damping  for  the  LMSR. 


26 


B35 


x  10*4  Heave/Pitch  Coupling  Added  Mass  -  LMSR 


Figure  5.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave/Pitch  Added  Mass  for  the 
LMSR. 


x  i  o"4  Heave/Pitch  Coupling  Damping  -  LMSR 


Figure  6.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave/Pitch  Damping  for  the  LMSR. 
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x  io"4  Pitch/Heave  Coupling  Added  Mass  -  LMSR 


Figure  7.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch/Heave  Added  Mass  for  the 
LMSR. 


x  io-4  Pitch/Heave  Coupling  Damping  -  LMSR 


Figure  8.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch/  Fleave  Damping  for  the 
LMSR. 


29 


x  1  o"4  Pitch  Added  Mass  -  LMSR 


Figure  1 1 .  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Added  Mass  for  the  Old 
LMSR  Data. 


Figure  12.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Damping  for  the  Old  LMSR 

Data. 
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Heave/Pitch  Coupling  Added  Mass-LMSR 


Figure  13.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave/Pitch  Added  Mass  for  the 


Old  LMSR  Data. 


Heave/Pitch  Coupling  Damping-LMSR 


Figure  14.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave/Pitch  Damping  for  the  Old 
LMSR  Data. 
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X104  Pitch/Heave  Coupling  Added  Mass  -  LMSR 


Figure  15.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch/Heave  Added  Mass  for  the  Old 

LMSR 


x  io"4  Pitch/Heave  Coupling  Damping  -  LMSR 


Figure  16.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave/Pitch  Damping  for  the  Old 
LMSR  Data. 
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x  •jo'fteave  Added  Mass  of  LMSR  due  to  Heave  Added  Mass  of  Barge 


Figure  21.  Plot  of  Multiple  Barge  Comparison  for 
Heave  Added  Mass  of  LMSR  due  to  Heave  Added 
Mass  of  Barge. 


x  lO-4  Heave  Damping  of  LMSR  due  to  Heave  Damping  of  Barge 


Figure  22.  Plot  of  Multiple  Barge  Comparison  for 
Heave  Damping  of  LMSR  due  to  Heave  Damping  of 

Barge. 
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x  i  o’7  Pitch  Added  Mass  of  LMSR  d  ue  to  Pitch  Added  Mass  of  Barge 


Figure  23.  Plot  of  Multiple  Barge  Comparison  for 
Pitch  Added  Mass  of  LMSR  due  to  Pitch  Added  Mass 
of  Barge. 

x  1 0  7  Pitch  Damping  of  LMSR  due  to  Pitch  Damping  of  Barge 


Figure  24.  Plot  of  Multiple  Barge  Comparison  for 
Pitch  Damping  of  LMSR  due  to  Pitch  Damping  of 
Barge. 
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Figure  25.  Plot  of  Multiple  Barge  Comparison  for 
Heave  Added  Mass  of  LMSR  due  to  Pitch  Added 
Mass  of  Barge. 


xIO*6  Heave  Damping  of  LMSR  due  to  Pitch  Damping  of  Barge 


Figure  26.  Plot  of  Multiple  Barge  Comparison  for 
Heave  Damping  of  LMSR  due  to  Pitch  Damping  of 
Barge. 
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x  i  o'^itch  Added  Mass  of  LMSR  due  to  Heave  Added  Mass  of  Barge 


Figure  27.  Plot  of  Multiple  Barge  Comparison  for 
Pitch  Added  Mass  of  LMSR  due  to  Heave  Added 
Mass  of  Barge. 


x  1  o"6  Pitch  Damping  of  LMSR  due  to  Heave  Damping  of  Barge 


Figure  28.  Plot  of  Multiple  Barge  Comparison  for 
Pitch  Damping  of  LMSR  due  to  Heave  Damping  of 
Barge. 
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Phase  Angle(in  Degrees)  Amplitude 


Plot  of  Wa\e  Exciting  Force  for  Heave  at  0  Degrees 


Figure  29.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  0  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  0  Degrees 


Figure  30.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle 
at  an  Incident  Wave  Angle  of  0  Degrees  for  the 
LMSR. 
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Plot  of  Wave  Exciting  Force  for  Pitch  at  0  Degrees 


Figure  3 1 .  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  0  Degrees  for  the  LMSR 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  0  Degrees 


Period(in  sec) 

Figure  32.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  0  Degrees  for  the  LMSR. 
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Plot  of  Wave  Exciting  Force  for  Heave  at  45  Degreess 


Figure  33.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  45  Degrees  for  the  LMSR 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  45  Degrees 


Figure  34.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  45  Degrees  for  the  LMSR. 
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Figure  35.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  45  Degrees  for  the  LMSR 


Plot  of  Ware  Exciting  Force  Phase  Angle  for  Pitch  at  45  Degrees 


Period(in  sec) 


Figure  36.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  45  Degrees  for  the  LMSR. 
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Phase  Angle(in  Degrees)  Amplitude 


Plot  of  Wave  Exciting  Force  for  Heave  at  90  Degrees 


Figure  37.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  90  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  90  Degrees 


Figure  38.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle 
at  an  Incident  Wave  Angle  of  90  Degrees  for  the 
LMSR. 
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x  io'3  Plot  of  Ware  Exciting  Force  for  Pitch  at  90  Degrees 


Figure  39.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  90  Degrees  for  the  LMSR. 


Plot  of  Ware  Exciting  Force  Phase  Angle  for  Pitch  at  90  Degrees 


Period(in  sec) 

Figure  40.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  90  Degrees  for  the  LMSR. 
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Phase  Angle(in  Degrees)  Amplitude 


Plot  of  Wa\£  Exciting  Force  for  Heave  at  135  Degrees 


Figure  41 .  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  135  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  135  Degrees 


Figure  42.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  135  Degrees  for  the  LMSR. 
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Figure  43.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  135  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  135  Degrees 


Period(in  sec) 

Figure  44.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  135  Degrees  for  the  LMSR. 
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Phase  Angle(in  Degrees)  Amplitude 


Plot  of  Wave  Exciting  Force  for  Heave  at  180  Degreess 


Figure  45.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  180  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  1 80  Degrees 


Figure  46.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  180  Degrees  for  the  LMSR. 
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Phase  Angle(in  Degrees)  ^  ^  Amplitude 


Plot  of  Wave  Exciting  Force  for  Pitch  at  1 80  Degrees 


Figure  47.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  180  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  1 80  Degrees 


Figure  48.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  1 80  Degrees  for  the 
LMSR. 
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Phase  Angle(in  Degrees)  Amplitude 


0.07 


Plot  of  Wave  Exciting  Force  for  Heave  at  225  Degrees 


Figure  49.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  225  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  225  Degrees 


Figure  50.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  225  Degrees  for  the 
LMSR. 
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Plot  of  Wave  Exciting  Force  for  Pitch  at  225  Degrees 


Figure  5 1 .  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  225  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  225  Degrees 


Figure  52.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of 225  Degrees  for  the 
LMSR. 
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Plot  of  Wave  Exciting  Force  for  Heave  at  270  Degreess 


Figure  53.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  270  Degrees  for  the  LMSR. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  270  Degrees 


Figure  54.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  270  Degrees  for  the 
LMSR. 
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Phase  Angle(in  Degrees) 


x  i  o'3  Plot  of  Wave  Exciting  Force  for  Pitch  at  270  Degrees 


Period(in  sec) 


Figure  55.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  270  Degrees  for  the 
LMSR. 

Plot  of  Wave  Exciting  Force  Phase  Angie  for  Pitch  at  270  Degrees 


0  2  4  6  8  10  12  14  16  18 

Period(in  sec) 

Figure  56.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  270  Degrees  for  the 
LMSR. 
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Plot  of  Wave  Exciting  Force  for  Heave  at  31 5  Degrees 


Figure  57.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  3 15  Degrees  for  the  LMSR 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  315  Degrees 


Figure  58.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Heave  Exciting  Force  Phase  Angle 
at  an  Incident  Wave  Angle  of  3 1 5  Degrees  for  the 
LMSR. 
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Figure  59.  Plot  of  Strip  Theory  and  Wamit  Program 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  3 1 5  Degrees  for  the  LMSR 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  31 5  Degrees 


Figure  60.  Plot  of  Strip  Theory  and  Wamit  Program  Comparisons  for  Pitch  Exciting 
Force  Phase  Angle  at  an  Incident  Wave  Angle  of  3 15  Degrees  for  the  LMSR. 
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B.  ROLL-ON,  ROLL-OFF  DISCHARGE  FACILITY  (BARGE)  ANALYSIS 

1.  Introduction 

The  hydrodynamic  characteristics  of  the  Roll-On,  Roll-Off  Discharge  Facility 
(RRDF)  are  presented  in  Figures  61  through  106.  These  figures  present  a  very  detailed 
analysis  of  the  barge  and  describe  the  effects  of  the  large  medium  speed  roll-on  roll-off 
vessel  (LMSR)  on  the  barge. 

2.  Data  Analysis 

Figure  61  illustrates  the  added  mass  in  heave  for  a  barge  unit  by  itself  and  also  the 
leading  (first)  and  trailing  (second)  barges  located  on  the  starboard  side  of  the  ship.  The 
results  displayed  in  all  the  figures  are  in  nondimensional  form  as  described  in  Table  1. 

It  should  be  noted  that  the  characteristic  length  which  is  utilized  to  nondimensionalize  the 
values  for  the  barges  are  that  of  the  ship  and  not  of  the  barge.  The  reasoning  for  this  is  to 
ensure  all  numbers  are  nondimensionalized  using  the  same  method,  which  allows  for  a 
more  accurate  comparison  of  results. 

The  results  displayed  are  those  for  the  three-dimensional  (3-D)  solver  only,  which 
modeled  the  RRDF  as  a  square  platform  with  shallow  draft.  The  actual  dimensions  for 
the  RRDF  are  given  in  Table  3.  Due  to  the  lack  of  a  precise  draft,  the  results  generated 
by  the  use  of  the  strip  theory  solver  were  determined  to  be  erroneous,  exhibiting  large 
divergences  from  the  3-D  solver  for  many  different  values.  Therefore,  data  was 
compared  against  a  stationary  barge  labeled  “1  Barge”  in  the  figures.  This  generated  the 
basis  for  establishing  the  relative  significance  of  hydrodynamic  proximity  effects.  The 
results  of  Figures  61  through  72  quantify  the  proximity  effects  on  the  heave  and  pitch 
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added  mass  and  damping.  These  figures  demonstrate  the  difference  from  the  baseline 
data  of  the  single  barge  to  the  barges  adjacent  to  the  LMSR  is  up  to  30  percent.  The 
irregular  oscillatory  behavior  of  the  plots  at  lower  periods  is  attributed  to  artificially 
generated  high  frequency  interactions  due  to  panel  discretization  effects. 

Results  from  the  wave  exciting  forces  are  shown  in  Figures  73  through  106. 
These  figures  display  the  non-dimensional  heaving  force  and  pitching  moment 
amplitudes  and  phase  angles,  for  the  headings  ranging  from  zero  to  360  degrees  at 
intervals  of  45  degrees.  Headings  of  zero  degrees  indicate  waves  coming  directly  from 
the  bow  of  the  RRDF,  90  degrees  from  the  starboard  beam,  180  degrees  from  the  stem, 
and  270  degrees  from  the  port  beam.  The  results  clearly  indicated  the  significance,  often 
two  or  more  times  higher,  of  the  proximity  affect  on  overall  exciting  forces.  The  very 
large  effect  of  proximity  forces  was  most  clearly  noted  in  the  following  figures: 

a.  Proximity  Induced  Exciting  Force 

Significant  effect  on  pitch  for  beam  seas  as  evidenced  in  Figure  83.  Here  the 
pitch  force  is  essentially  zero  for  the  single  barge  case  which  is  expected  due  to  hull 
geometry.  Therefore,  any  non-zero  forces  acting  on  the  barges  are  entirely  due  to  the 
presence  of  the  LMSR.  Similar  results  are  found  when  the  wave  angle  is  from  270 
degrees  in  Figure  99. 

b.  Sheltering/Magnification  Effect 

It  should  be  noted  that  both  the  heave  and  pitch  exciting  forces  are  significantly 
larger  when  the  waves  are  from  90  degrees  compared  to  when  the  waves  are  from  270 
degrees.  The  reason  for  this  disparity  is  due  to  the  incident  waves  being  reflected  off  the 
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starboard  side  of  the  LMSR  and  then  impacting  upon  the  RRDF’s  which  again  create  a 
magnification  effect  in  the  90  degree  case.  By  contrast,  when  the  wave  direction  is  from 
270  degrees  the  incident  waves  are  scattered  by  the  port  side  of  the  ship  that  produces  a 
sheltering  effect  on  the  RRDF’s  from  wave  action. 
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3.  Figures 


Period(in  seconds) 


Figure  61 .  Plot  of  Single  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Added  Mass. 
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Period(in  seconds) 

Figure  62.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Damping. 
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Figure  63.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Added  Mass. 


Period(in  seconds) 


Figure  64.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Damping. 
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x  i  o'7  Pitch  Added  Mass  of  Barge  due  to  Pitch  Added  Mass  of  LMSR 


Figure  67.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Added  Mass  of  Barge  Due  to  Pitch  Added  Mass 
of  LMSR. 


x  io~  Pitch  Damping  of  Barge  due  to  Pitch  Damping  of  LMSR 
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Figure  68.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Damping  of  Barge  Due  to  Pitch  Damping  of 
LMSR. 
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Period(in  seconds) 


x  •)  o"3  Plot  of  Wave  Exciting  Force  for  Heave  at  0  Degrees 


Figure  73.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  0  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  0  Degrees 


Figure  74.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  Phase  Angle 
at  an  Incident  Wave  Angle  of  0  Degrees. 
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x  -jo  5  Plot  of  Wave  Exciting  Force  for  Pitch  at  0  Degrees 


Figure  75.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  0  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  0  Degrees 


Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  0  Degrees. 
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Period(in  sec) 


Figure  77.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  45  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  45  Degrees 


Period(in  sec) 


Figure  78.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  45  Degrees. 


x  i  o"5  Plot  of  Wave  Exciting  Force  for  Pitch  at  45  Degrees 


Figure  79.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  45  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  45  Degrees 


Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  45  Degrees. 
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Period(in  sec) 


Figure  81.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  90  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  90  Degrees 
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Figure  82.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  90  Degrees. 
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x  10*4  Plot  of  Wave  Exciting  Force  for  Pitch  at  90  Degrees 


Figure  83.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  90  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  90  Degrees 


Period(in  sec) 

Figure  84.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  90  Degrees. 
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x  i  o’  Plot  of  Wave  Exciting  Force  for  Heave  at  1 35  Degreess 


Figure  85.  Plot  of  Multiple  Barge  Comparisons  for 
Heave  Exciting  Force  at  an  Incident  Wave  Angle  of 
135  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  135  Degrees 


Figure  86.  Plot  of  Multiple  Barge  Comparisons  for 
Heave  Exciting  Force  Phase  Angle  at  an  Incident 
Wave  Angle  of  135  Degrees. 
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Plot  of  Wave  Exciting  Force  for  Pitch  at  135  Degrees 


Figure  87.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Exciting  Force  at  an  Incident  Wave  Angle  of  135 

Degrees. 

Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  135  Degrees 


Figure  88.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Exciting  Force  Phase  Angle  at  an  Incident  Wave 
Angle  of  135  Degrees. 
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Plot  of  Wave  Exciting  Force  for  Heave  at  180  Degrees 
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Figure  89.  Plot  of  Multiple  Barge  Comparisons  for 
Heave  Exciting  Force  at  an  Incident  Wave  Angle  of 
180  Degrees. 

Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  180  Degrees 
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Figure  90.  Plot  of  Multiple  Barge  Comparisons  for 
Heave  Exciting  Force  Phase  Angle  at  an  Incident 
Wave  Angle  of  180  Degrees. 
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x  1 0  5  Plot  of  Wave  Exciting  Force  for  Pitch  at  1 80  Degrees 


Figure  91.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Exciting  Force  at  an  Incident  Wave  Angle  of  1 80 

Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  180  Degrees 


Figure  92.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Exciting  Force  Phase  Angle  at  an  Incident  Wave 
Angle  of  180  Degrees. 
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Plot  of  Wave  Exciting  Force  for  Pitch  at  225  Degrees 


Figure  95.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Exciting  Force  at  an  Incident  Wave  Angle  of 
225  Degrees. 

Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  225  Degrees 


Figure  96.  Plot  of  Multiple  Barge  Comparisons  for 
Pitch  Exciting  Force  Phase  Angle  at  an  Incident 
Wave  Angle  of  225  Degrees. 
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x  i  o  Plot  of  Wave  Exciting  Force  for  Heave  at  270  Degrees 


Figure  97.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  270  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  270  Degrees 


Figure  98.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  Phase  Angle 
at  an  Incident  Wave  Angle  of  270  Degrees. 
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Plot  of  Wave  Exciting  Force  for  Pitch  at  270  Degrees 


Figure  99.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  270  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  270  Degrees 


Figure  100.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  270  Degrees. 
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x  •)  o’3  Plot  of  Wave  Exciting  Force  for  Heave  at  31 5  Degrees 


Figure  101.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  315  Degrees. 


Plot  of  Wave  Exciting  Force  Phase  Angle  for  Heave  at  315  Degrees 


Period(in  sec) 


Figure  102.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  3 15  Degrees. 
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Plot  of  Wave  Exciting  Force  for  Pitch  at  315  Degrees 


Figure  103.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  3 1 5  Degrees. 

Plot  of  Wave  Exciting  Force  Phase  Angle  for  Pitch  at  315  Degrees 


Figure  104.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  Phase  Angle  at 
an  Incident  Wave  Angle  of  315  Degrees. 
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x  io'3Plot  of  Wave  Exciting  Force  for  Heave  at  90  and  270  Degreess 
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Figure  105.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Heave  Exciting  Force  at  an  Incident 
Wave  Angle  of  90  Degrees  and  270  Degrees. 


x  i  o"4  Plot  of  Wave  Exciting  Force  for  Pitch  at  90  and  270  Degrees 


Figure  106.  Plot  of  One  Barge  and  Multiple  Barge 
Comparisons  for  Pitch  Exciting  Force  at  an  Incident 
Wave  Angle  of  90  Degrees  and  270  Degrees. 
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C.  RAMP  ANALYSIS 


1.  Introduction 

The  forces  on  the  ramp  exerted  from  the  hydrodynamic  motions  of  the  LMSR  and 
barge  are  presented  in  Figures  1 07  through  131a.  These  figures  show  the  overall  forces 
on  the  pinned  locations  that  connect  the  ramp  to  the  LMSR  and  the  barge.  To  derive 
these  forces,  the  data  utilized  was  that  obtained  from  WAMIT  as  this  data  produced 
excellent  results  when  compared  with  actual  wave  tank  model  tests.  The  ramp  is 
assumed  to  be  a  rigid  body  for  analysis  and  stiffness  of  the  spring  connections  is  assumed 
to  be  proportional  to  the  hydrostatic  heave  stiffness  of  the  LMSR  or  barge,  depending  on 
which  pinned  connection  is  being  analyzed.  Table  4  shows  the  dimensions  of  the  ramp. 


Description 

Dimensions 

Ramp  Length 
(in  Feet,  Pivot  to  Pivot) 

165 

Ramp  Width  (in  Feet,  Overall) 

20 

Ramp  Weight 

96.94  Long  Tons 

Table  4.  Ramp  Dimensions 


2.  Data  Analysis 

Figures  107  through  122  show  a  plot  of  the  ramp  force  magnitudes  for  the  three- 
pinned  connections.  It  is  clear  that  the  pinned  connection  that  receives  the  largest  peak 
magnitude  force  is  the  third  king  pin  which  connects  the  ramp  to  the  barge.  The 
connection  that  receives  the  largest  force  throughout  the  entire  range  is  the  first  pinned 
connection.  The  third  connection  receives  up  to  5  times  the  magnitude  of  force  as 
compared  to  the  pinned  connections  attaching  the  ramp  to  the  LMSR  at  the  peak 
frequency.  The  first  king  pin  has  a  force  about  1 .5  times  greater  than  the  third  connection 


82 


throughout  the  entire  range.  The  forces  on  the  barge’s  king  pin  are  largest  when  the 
waves  are  from  90  and  270  degrees. 

When  looking  at  the  individual  king  pin  for  the  LMSR  and  the  barge, 
comparisons  were  made  between  the  actual  force  at  the  pinned  connection  and  the 
hydrodynamic  heave  force  acting  on  the  vessel.  This  data  is  displayed  in  Figures  123, 
125,  and  127  where  the  actual  force  which  each  king  pin  experiences  is  higher  than  the 
hydrodynamic  heave  force  on  the  vessel  at  higher  wave  frequencies.  The  force  that  the 
ramp  connections  experience  drops  off  at  longer  periods. 

Figures  107  through  122  and  129  through  131a  show  the  changes  in  overall 
excitation  force  on  the  pinned  connection  if  different  spring  stiffnesses  are  assumed  to 
exist.  As  can  be  seen  in  Figures  129, 130,  and  131,  the  larger  the  value  of  the  spring 
constant,  the  greater  the  magnitude  of  force  which  each  king  pin  observes.  This  implies 
that  a  lower  stiffness  value  for  the  spring  connection  is  one  method  of  minimizing 
excitation  forces  on  the  pinned  connections. 
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Magnitude  of  Force 


3.  Figures 

Plot  of  Force  Magnitudes  at  Various  Angles  for  1st  Pinned  Connection,  K=0.01C33 


Figure  107.  Plot  of  Ramp  Exciting  Force  on  the  First 
Pinned  Connections  at  Various  Wave  Angles  with  a 
Spring  Constant  of  0.01 . 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  1st  Pinned  Connection 


Figure  108.  Plot  of  Ramp  Exciting  Force  on  the  1st 
Pinned  Connection  at  Various  Wave  Angles  with  a 
Spring  Constant  of  0.1. 
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Magnitude  of  Force  ^  Magnitude  of  Force 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  1st  Pinned  Connection,  K=C33 
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Figure  109.  Plot  of  Ramp  Exciting  Force  on  the  First 
Pinned  Connections  at  Various  Wave  Angles  with  a 
Spring  Constant  of  1. 


Figure  110.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connections  at  Various  Wave  Angles 
with  a  Spring  Constant  of  0.0 1 . 
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Magnitude  of  Force  Magnitude  of  Force 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  2nd  Pinned  Connection,  K=0.01C33 


Figure  1 10a.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connections  at  Various  Wave  Angles 
with  a  Spring  Constant  of  0.01. 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  2nd  Pinned  Connection 


Figure  111.  Plot  of  Ramp  Exciting  Force  on  the  2nd 
Pinned  Connection  at  Various  Wave  Angles  with  a 
Spring  Constant  of  0. 1 . 
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Magnitude  of  Force 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  2nd  Pinned  Connection,  K=  1C33 


Figure  111a.  Plot  of  Ramp  Exciting  Force  on  the  2nd 
Pinned  Connection  at  Various  Wave  Angles  with  a 
Spring  Constant  of  0. 1 . 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  2nd  Pinned  Connection,  K=C33 


Figure  1 12.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connection  at  Various  Wave  Angles 
with  a  Spring  Constant  of  1. 
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Magnitude  of  Force 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  2nd  Pinned  Connection,  K=C33 


Figure  1 12a.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connection  at  Various  Wave  Angles 
with  a  Spring  Constant  of  1 . 

Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  3rd  Pinned  Connection,  K=0.01C33 


Figure  113.  Plot  of  Ramp  Exciting  Force  on  the 
Third  Pinned  Connection  at  Various  Wave  Angles 
with  a  Spring  Constant  of  0.0 1 . 
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Magnitude  of  Force 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  3rd  Pinned  Connection,  K=0.01C33 


Figure  113a.  Plot  of  Ramp  Exciting  Force  on  the 
Third  Pinned  Connection  at  Various  Wave  Angles 
with  a  Spring  Constant  of  0.01. 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  3rd  Pinned  Connection 


Figure  1 14.  Plot  of  Ramp  Exciting  Force  on  the  3rd 
Pinned  Connection  at  Various  Wave  Angles  with  a 
Spring  Constant  of  0.1 . 
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Magnitude  of  Force 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  3rd  Pinned  Connection,  K=  1C33 


Figure  1 14a.  Plot  of  Ramp  Exciting  Force  on  the  3rd 
Pinned  Connection  at  Various  Wave  Angles  with  a 
Spring  Constant  of  0.1 . 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  3rd  Pinned  Connection,  K=C33 


Figure  115.  Plot  of  Ramp  Exciting  Force  on  the  Third 
Pinned  Connections  at  Various  Wave  Angles  with  a 
Spring  Constant  of  1 . 
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Magnitude  of  Force  Magnitude  of  Force 


Plot  of  Ramp  Force  Magnitudes  at  Various  Angles  for  3rd  Pinned  Connection,  K=C33 
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Figure  116.  Plot  of  Ramp  Exciting  Force  on  all  Three 
Pinned  Connections  at  a  Wave  Angle  of  90  Degrees 
with  a  Spring  Constant  of  0. 1 . 


Figure  1 1 5a.  Plot  of  Ramp  Exciting  Force  on  the 
Third  Pinned  Connections  at  Various  Wave  Angles 
with  a  Spring  Constant  of  1 . 


Plot  of  Force  on  the  ramp  at  90  degrees 
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Force  on  all  Spring  Connections  at  90  degrees,  K  that  of  0.1  ofC33 


Figure  116a.  Plot  of  Ramp  Exciting  Force  on  all 
Three  Pinned  Connections  at  a  Wave  Angle  of  90 
Degrees  with  a  Spring  Constant  of  0. 1 . 


Plot  of  Phase  on  the  ramp  at  90  degrees 


Figure  117.  Plot  of  Ramp  Exciting  Force  Phase 
Angle  on  all  Three  Pinned  Connections  at  a  Wave 
Angle  of  90  Degrees  with  a  Spring  Constant  of  0. 1 . 
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Magnitude  of  Force 


0.25 


Force  on  all  Spring  Connections  at  90  degrees,  K  that  of  0.1  of  LMSR  C33 


Figure  118.  Plot  of  Ramp  Exciting  Force  on  all 
Pinned  Connections  at  a  Wave  Angle  of  90  Degrees 
with  all  Spring  Constants  0.1  of  the  LMSR  Stiffness. 


Force  on  all  Spring  Connections  at  90  degrees,  K  that  of  0.1  ofLMSRC33 


Figure  118a.  Plot  of  Ramp  Exciting  Force  on  all 
Pinned  Connections  at  a  Wave  Angle  of  90  Degrees 
with  all  Spring  Constants  0.1  of  the  LMSR  Stiffness. 
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Plot  of  Force  on  the  ramp  at  1 35  degrees 


Figure  119.  Plot  of  Ramp  Exciting  Force  on  all  Three 
Pinned  Connections  at  a  Wave  Angle  of  135  Degrees 
with  a  Spring  Constant  of  0. 1 . 


Plot  of  Force  on  the  ramp  at  1 80  degrees 


Figure  120.  Plot  of  Ramp  Exciting  Force  on  all  Three 
Pinned  Connections  at  a  Wave  Angle  of  180  Degrees 
with  a  Spring  Constant  of  0. 1 . 
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Magnitude  of  Force  Magnitude  of  Force 


Plot  of  Force  on  the  ramp  at  225  degrees 


Figure  121.  Plot  of  Ramp  Exciting  Force  on  all  Three 
Pinned  Connections  at  a  Wave  Angle  of  225  Degrees 
with  a  Spring  Constant  of  0. 1 . 


Plot  of  Force  on  the  ramp  at  270  degrees 


Figure  122.  Plot  of  Ramp  Exciting  Force  on  all  Three 
Pinned  Connections  at  a  Wave  Angle  of 270  Degrees 
with  a  Spring  Constant  of  0.1. 


95 


Phase  Angle  Magnitude  of  Force 


Plot  of  LMSR  Ramp  Force  Magnitude  and  LMSR  Heave  Motion  at  90  Degrees 


Figure  123.  Plot  of  Ramp  Exciting  Force  on  the  First 
Pinned  Connection  and  Ship  Heave  Motion  at  a  Wave 
Angle  of  90  Degrees  with  a  Spring  Constant  of  0.1. 


Plot  of  LMSR  Ramp  Force  Phase  Angle  and  LMSR  Heave  Phase  Angle  at  90  Degrees 


Figure  124.  Plot  of  Ramp  Exciting  Force  on  the  First 
Pinned  Connection  and  Ship  Heave  Motion  Phase 
Angles  at  a  Wave  Angle  of  90  Degrees  with  a  Spring 
Constant  of  0.1. 
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Plot  of  Barge  Ramp  Force  Magnitude  and  Barge  Heave  Motion  at  90  Degrees 


Figure  125.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connection  and  Ship  Heave  Motion  at 
a  Wave  Angle  of  90  Degrees  with  a  Spring  Constant 

of  0.1. 

Plot  of  Barge  Ramp  Force  Phase  Angle  and  Barge  Heave  Phase  Angle  at  90  Degrees 


Figure  126.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connection  and  Ship  Heave  Motion 
Phase  Angles  at  a  Wave  Angle  of  90  Degrees  with  a 
Spring  Constant  of  0. 1 . 
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Phase  Angle  Magnitude  of  Force 


Plot  of  Barge  Ramp  Force  Magnitude  and  Barge  Heave  Motion  at  90  Degrees 


Figure  127.  Plot  of  Ramp  Exciting  Force  on  the  Third 
Pinned  Connection  and  Ship  Heave  Motion  at  a  Wave 
Angle  of  90  Degrees  with  a  Spring  Constant  of  0.1 . 

Plot  of  Ramp  Force  Phase  Angle  and  Barge  Heave  Phase  Angle  at  90  Degrees 


Figure  128.  Plot  of  Ramp  Exciting  Force  on  the  Third 
Pinned  Connection  and  Ship  Heave  Motion  Phase 
Angles  at  a  Wave  Angle  of  90  Degrees  with  a  Spring 
Constant  of  0.1. 
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Plot  of  Force  on  the  ramp  at  90  degrees,  Various  K  Values 


Figure  129.  Plot  of  Ramp  Exciting  Force  on  the  First 
Pinned  Connection  at  a  Wave  Angle  of  90  Degrees 
with  Various  Spring  Constants. 


Plot  of  Force  on  the  ramp  at  90  degrees,  2nd  Pinned  Connection,  Various  K  Values 
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Figure  130.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connection  at  a  Wave  Angle  of  90 
Degrees  with  Various  Spring  Constants. 
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Magnitude  of  Force 


Plot  of  Force  on  the  ramp  at  90  degrees,  2nd  Pinned  Connection,  Various  K  Values 


Figure  130a.  Plot  of  Ramp  Exciting  Force  on  the 
Second  Pinned  Connection  at  a  Wave  Angle  of  90 
Degrees  with  Various  Spring  Constants. 


Plot  of  Force  on  the  ramp  at  90  degrees,  3rd  Pinned  Connection,  Various  K  Values 


Figure  1 3 1 .  Plot  of  Ramp  Exciting  Force  on  the 
Third  Pinned  Connections  at  a  Wave  Angle  of  90 
Degrees  with  Various  Spring  Constants. 
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V.  CONCLUSIONS 


The  main  conclusions  from  this  work  are  summarized  below: 

1 .  The  comparison  between  a  panel  based  three-dimensional  hydrodynamic 
solver  and  a  two-dimensional  strip  theory  solver  showed  excellent  agreement  when 
modeling  the  LMSR.  When  the  LMSR  was  modeled  with  the  barges  attached,  although 
the  strip  theory  solver  was  unable  to  account  for  the  wave  interaction  affects,  the  added 
mass  and  damping  coefficients  as  well  as  the  wave  exciting  forces  for  the  LMSR  were 
still  predicted  fairly  well. 

2.  The  use  of  the  strip  theory  solver  is  more  economical  in  that  it  requires 
significantly  less  computational  time.  Therefore,  it  may  be  effectively  used  to  model  the 
LMSR.  A  combination  of  two-dimensional  and  three-dimensional  hydrodynamic  solvers 
is  the  most  efficient  method  to  obtain  the  added  mass  and  damping  values  for  the 
LMSR/barge  system. 

3.  The  strip  theory  solver  was  unable  to  accurately  predict  the  hydrodynamic 
effects  of  the  LMSR  on  the  barge.  This  was  attributed  mainly  to  the  irregular 
frequencies. 

4.  Exciting  forces  that  are  induced  due  to  proximity  produce  a  pitching 
moment  on  the  barges  attached  to  the  LMSR  in  beam  seas.  This  moment  causes  the 
barges  to  pitch  when  they  would  naturally  not  exhibit  any  pitch  motion.  This  results  in 
larger  barge  motion  than  would  be  expected  in  beam  seas  and  thus  produces  a  large  force 
on  the  pinned  connections  in  beam  seas. 
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5.  A  sheltering  and  magnification  effect  produces  a  dramatically  lowered  total 
force  on  the  barges  when  seas  are  from  270  degrees  and  a  significantly  larger  force  when 
the  wave  direction  is  from  90  degrees. 

6.  The  king  pin  on  the  LMSR  is  the  point  where  the  most  force  is  exhibited. 
This  point  is  anywhere  from  2  to  3  times  larger  than  the  connection  on  the  barge.  At 
higher  frequencies,  this  force  may  be  higher  than  the  hydrodynamic  forces  which  the  ship 

experiences. 

7.  With  larger  spring  constants,  the  force  acting  on  the  pinned  connections  is 
significantly  greater. 


104 


VI.  RECOMMENDATIONS 


Based  upon  the  data  and  conclusions,  the  following  recommendations  are  made: 

1 .  Reformulate  the  problem  as  a  four  degree  of  freedom  system.  This  would 
require  modeling  the  ramp  as  non-rigid  body.  Preliminary  studies  indicate  that  the  first 
torsional  mode  of  the  ramp  is  the  most  significant,  so  this  would  need  to  be  incorporated 
first  in  the  analysis. 

2.  Utilize  the  transfer  functions  that  were  developed  in  this  work  to  select  the 
proper  isolator  properties.  The  minimization  index  would  be  the  maximum  stress  in  the 
ramp. 

3.  Analyze  all  six  motions  of  the  LMSR  and  barges  to  verify  that  other  motions 
are  negligible  as  assumed  in  this  problem. 
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APPENDIX  A 


FIGURE  OF  LMSR/BARGE  CONNECTION: 


Figure  A.  Basic  Drawing  of  LMSR  and  Roll-On,  Roll-Off  Discharge  Facility  Connected 
together  by  Ramp.  Black  Dots  Indicate  the  Locations  and  Number  of  the  Pinned 

Connections. 
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